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^ . ABSTRACT 
^ . 

^ . We have investigated the structure and evolution of hot gas in the 30 Dor nebula, 

O ■ based on recent X-ray observations. Our deep ROSAT lllil image shows that diffuse 

! X-ray emission arises in blister-shaped regions outlined by loops of HII gas. X-ray 

CN| I spectroscopic data from ASCA confirm the thermal nature of the emission and indicate 

that hot gas temperature decreases from the core to the halo of the nebula. The 
structure of the nebula can be understood as outflows of hot and HII gases from the 
parent giant molecular cloud of the central OB association. The dynamic mixing 
C3 ' between the two gas phases is likely responsible for the mass loading to the hot gas, as 

required to explain the observed thermal structure and X-ray luminosity of the nebula. 
OO '. Such processes should also be important in the formation of similar giant HII regions 

and in their subsequent evolution into supergiant bubbles or galactic chimneys. 



O ■ Subject headings: ISM: individual (30 Dor) — ISM: bubbles — ISM: structure — ISM: 

> ' HII regions — Magellanic Clouds — galaxies: ISM — X-ray: ISM 



X 



1. Introduction 

As the most luminous HII nebula in the Local Group of galaxies, 30 Dor in the Large 
Magellanic Cloud (LMC) provides the most accessible paradigm to understand massive star 
formation and its interaction with the interstellar medium (ISM). Responsible for the ionization 
of this spectacular nebula 10^^ Ly photon s^^) is the OB association NGC 2070. The central 
cluster R136 (r ^ 10 pc; assuming the LMC distance as D=47 kpc; Gould 1995) alone accounts for 
about half of the required ionization radiation (Walborn & Blades 1997; Crowther & Dessart 1998 
and references therein). The bulk of very massive stars, contained primarily within the compact 
core of the cluster (r ~ 2 pc), were formed in a star forming burst about 1-2 Myrs ago (Massey &: 
Hunter 1998). HII filaments form shells or loops that often extend more than 100 pc away from 
the central cluster (e.g., Chu & Kennicutt 1994). However, no model has yet been proposed to 
explain the origin and evolution of this nebula structure. 

The 30 Dor nebula is the first giant HII region (GHR) shown to be a strong diffuse X-ray 
emitter. As seen in an IPC image of the Einstein Observatory (Wang & Helfand 1991), the 
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emission peaks within HII shells, indicating the presence of hot gas of a few times 10^ K. Stellar 

winds from massive stars and supernova explosions are presumably responsible for the heating of 
the gas. Relatively recent observations with the ROSAT Positional Sensitive Proportional Counter 
(PSPC) have confirmed the Einstein results (Norci & Ogelman 1993; Chu 1993). 

In this Letter, we present a study of 30 Dor based chiefly on a deep X-ray image from the 
ROSAT High Resolution Imager (RHRI) and on an X-ray spectrum from the Advanced Satellite 
for Cosmology and Astrophysics {ASCA). These two complimentary data sets enable us to 
spatially resolve the detailed X-ray emission structure and to characterize the thermal properties 
of the hot gas. We compare the X-ray data with optical and UV observations and extend our 
discussion (Wang 1996) about the origin of the hot gas and its subsequent evolution. 

2. Description of the X-ray Data 

The RHRI image, as presented in Fig. 1, is a co-add of three observations: ROSAT ^ 
rh600228n00 (30 ks exposure), rh400779a01 (79 ks), and rh400779n00 (26 ks). The first observation 
has already been used in the study of the two bright point-like sources in the core of the nebula, 
which are Wolf-Rayet -I- black hole binary candidates (Fig. 2; R140 and Mrk 34; Wang 1995), and 
the Crab- like supernova remnant N157B (Wang & Gotthelf 1998a). The second observation has 
been utilized for both timing and positioning of the recently discovered 16 ms pulsar within the 
remnant (Wang & Gotthelf 1998b and references therein). Briefly, the RHRI data have a spatial 
resolution ranging from ~ 6" (FWHM) near the image center to 30" near the edges, or a factor 
of up to ~ 10 (5) better than the resolution of the IPC (PSPC). The energy coverage is in the 
range of 0.1 - 2 keV. We processed the RHRI data with a software provided by S. Snowden for 
data editing, background subtraction, and exposure correction. We then adaptively smoothed 
the resultant X-ray intensity image with a Gaussian of adjustable size to achieve a uniform 
signal-to-noise ratio of 6 for Fig. 1 and 4 for Fig. 2. 

The X-ray spectral data were obtained from the ASCA SIS observation # ad20000000. Wang 
& Gotthelf (1998a) have presented two broad band images constructed from this observation, which 
had a spatial resolution of ~ 1' (FWHM) and a spectral resolution of 6E/E ~ 0.02(5.9 keV/E)^-^. 
The observation, taken in a 4-CCD mode, was pointed at R.A., Dec. (J2000)= 5'*38™32^-69°9'44" 
and had a roll angle of 20°. 3 north to east. We utilized only the data from the two northern 
CCD chips of each sensor to minimize the contamination caused by scattered counts from 
N157B, PSR B0540-69, and LMC X-1 to the south. Specifically, the on-nebula spectrum was 
extracted from a half circle of 7^5 radius, centered at the pointing direction, plus a northern 
stretch of 5^4 wide to complete the coverage of the northern X-ray spur of the nebula (Fig. 1). 
We estimated the background from the same on-nebula chip areas in the adjacent observation 
(ad90001000 center at 5''43'"31^-69°13'48"), just east of 30 Dor. We also tested an alternative 
background estimate from the regions away from the nebula, but in the same on-nebula CCD 
chips. The resultant background-subtracted spectra are essentially the same; the diflFerences 
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in best-fit spectral parameters (Table 1) are all within 5%. We added the spectra from both 
sensors of each observation (on- or off-nebula) for subsequent spectral model fitting. This spectral 
data reduction followed the procedure prescribed in the on-line ASCA Data Reduction Guide 
([http://h easarc.gsfc.nasa.gov /docs/ asca/ abc /abc.html| ) . 

3. Results 

The X-ray emission from 30 Dor is predominantly diffuse in origin. Fig. 1 shows no evidence 
for point-like X-ray sources which may be related to 30 Dor, except for the two X-ray binary 
candidates mentioned above (§2) and the apparent X-ray contribution from R136 (Fig. 2). In 
general, there is no detailed correlation between the X-ray emission and the presence of massive 
stars, as appear in the Astro-1 UIT UV image (Fig. 3; Hill et al. 1993). The detection limit 
for an individual point-like source is ~ 1 x 10^^ ergs s~^ (0.5-2 keV). Enhanced diffuse X-ray 
emission is enclosed within HII loops, which are typically anchored to dense molecular clouds in 
the core of the nebula (Wang 1996; Fig. 3). But the detailed morphology of individual features 
remains uncertain. It is still difficult to disentangle various projection effects and differential 
X-ray absorption across the field. For example, the two shell-like features of ~ 3' diameters to the 
south of R136 may well be the projection of two loops that originate in the core of the nebula. 
Apparently, the 30 Dor nebula is a complex of diffuse X-ray-emitting blisters. 

In the central region (Fig. 2), strong X-ray emission arises in cavities outlined by ionization 
fronts (Scowen et al. 1998) at boundaries of molecular clouds (e.g., Wang 1996; Johansson et al. 
1998; Rubio et al. 1998). The cavity around R136, in particular, is a well-defined blister, which 
is open, at least, to the east (e.g., Dickel et al. 1994). HII gas kinematics (e.g, Chu &: Kennicutt 
1994) show that the eastern part of this blister is expanding much more rapidly than the west 
part. This expansion is most likely driven by X-ray-emitting gas that is flowing out from near 
R136 to the halo of the nebula. The transition zone between the core and the halo of the nebula 
appears at r ~ 2' — 4' from R136, at outer boundaries of the parent giant molecular cloud (GMC) 
of the OB association. Particularly in the east-west direction, Ha-emitting filaments clearly show 
the morphology of outward streamers (Fig. 3). 

The ASCA spectral data (Fig. 4) confirm that X-ray emission in the 0.5-2 keV band is 
predominantly thermal in origin. Emission lines are clearly present: e.g., Ne 1.0 keV), Mg 
(1.3 keV), and Si (1.9 keV). The flat spectrum at higher energies, however, suggests the presence 
of a nonthermal component, which may represents the expected hard X-ray spectral tails of the 
black hole binaries (Wang 1995). A power law with a typical photon index of 2 describes this 
component well, although its exact spectral shape is uncertain due to the limited energy coverage 
and statistics of the data. Using the Raymond &: Smith optically thin plasma model to characterize 
the thermal portion of the ASCA spectrum, we find that at least two temperature components are 
needed for a reasonable good fit to the overall spectral shape. But, this three-component model 
is still an oversimplification. Even when the abundances are allowed to adjust, the model does 
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not adequately account for several line features which peak at intermediate plasma temperatures 
of a few times 10^ K. The inclusion of these lines in the model (Table 1) gives an acceptable fit 
to the data {x^ /n.d.f. = 124/107; Fig. 4). The two thermal components characterize a minimum 
temperature range of the X-ray-emitting gas, ~ 2 — 9 x 10® K. The high temperature component 
dominates in the energy band 1.5 keV, and arises primarily in the core region around R136 
(r ;^ 1'; Wang & Gotthelf 1998a). Howvever, this component cannot be explained by the central 
two-point hke sources (Fig. 2); their total RHRI count rate is 0.013 ± 0.001 (Wang 1995), or only 
a quarter of the predicted count rate of the component alone. Therefore, the hard temperature 
component also originates predominately in diffuse hot gas and its temperature apparently 
decreases from the core to the halo of the nebula. 

The X-ray spectral analysis further provides a measurement of the mean X-ray-absorbing 
gas column density toward 30 Dor. A comparison of the measured column density Nh (Table 
1) with the mean reddening E(5 - V) = 0.4 - 0.5 of 30 Dor (e.g., Parker 1993; Dickel et 
al. 1994) gives Nh/E(B — F) 3.1 x 10^^ cm~^ mag~^, which is considerably greater than 
Nhi/E(B — y) « 2.4 X 10^^ cm~^ mag~^ (Fitzpatrick 1986). This is expected because Nh includes 
contributions from both partially-ionized and molecular gas phases, in addition to the H I column 
density Nhi- The X-ray absorption in the ASCA energy range (^0.5 keV) is insensitive to the 
ionization and chemical states of the gas (Morrison & McCammon 1983). 

Using the best-fit spectral model (Table 1), we estimate physical parameters of hot gas in the 
nebula. An approximate conversion between the observed RHRI intensity to the emission measure 
of hot gas is ~ 40 cvoT^ pc/(counts ks"-*^ arcmin"^) with an uncertainty of ~ 2 within the ranges 
of the spectral parameters. The total count rate from a circle of 10' radius around R136, excluding 
the N157B and point-like source contributions, is ~ 0.26 counts s~^. This rate corresponds to 
an intrinsic luminosity La;(0-5 — 2 keV) ~ 9 x 10^^ ergs s~^ and a radiative cooling rate of 
~ 1 X 10^^ ergs s~^. The spectral parameters of the high and low temperature thermal components 
suggest an average pressure p/k ~ 2 X 10^ K cm ^ in the core and a factor of ~ 7 lower in the 
halo of the nebula. The total thermal energy and mass of the hot gas, if within a sphere of 
r = 140 pc and with a volume filling factor /0.2 (in units of 20%), are ~ (3 X 10^2 grgs) ji/2 
(4 X 10^Mo)/o 2^, contained mostly in the low temperature component. 



4. Discussions 

The formation of the giant 30 Dor nebula is clearly driven by the intense energy release 
from massive stars. With a total mass losing rate of ~ 1O~^M0 yr~^ and a typical terminal 
wind velocity of ~ 3 x 10'^ km s~^, the central cluster R136 alone has a stellar wind luminosity 
of ~ 3 X 10^^ ergs s~^ (e.g., Chu &: Kennicutt 1994). Over its lifetime, the cluster should have 
released ~ 10^^ ergs of mechanical energy. A substantial fraction of this energy is contained in 
the hot gas, and that the radiative cooling of the hot gas is considerable (§3). The remaining 
energy could be accounted for by the global and turbulent motion of various gas components (e.g.. 
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Chu & Kennicutt 1994) and by the energy loss to the formation of dense gas shells or filaments 
(Mac Low & McCray 1988). In comparison, the hot gas in the central cavity accounts for only 
~ 10^^ ergs. Thus, outflows must have occurred. Through a highly inhomogeneous and clumpy 
medium, the outflows can naturally lead to the formation of X-ray-emitting blisters of various 
shapes and kinematics (Wang 1996). 

The key process invloved in the evolution of the hot gas is mass loading. Just behind the 
terminal shock, stellar wind materials have a temperature of ~ 10^ K, which is one order of 
magnitude greater than the measured average temperature of hot gas even in the core of the 
nebula. The mass loading can naturally increase the density and decrease the temperature of the 
materials, neccessary for explaining the observed X-ray emission. 

The most effective mass loading process is likely the dynamic mixing of hot gas with 
HII gas. While the mass efficiency of massive star formation is typically only 5-10% (e.g., 
McKee 1989), much of the parent GMC of the OB association is being eroded via photon- 
evaporation. The erosion has been accellerated greatly since the nebula broke out from 
the GMC. The mass eroding rate can be estimated as (Whitworth 1979; York et al. 1989) 
~ O.12M0 yr-i(5i/10^2 Ly photons s-^)'^^^U/W^ yif^'^ {uo/lO^ cm-^)-^^^ where Si, U, and 
are the ionizing flux of NGC 2070, the effective age of the OB association with the present flux, 
and the mean density of the GMC. Only if ^ 10% of this photon-evaporated HII gas is loaded 
to the shocked wind materials in the core of the nebula, can the temperature of hot gas there be 
explained. Over the lifetime of the nebula, the total evaporated mass is ~ 10^ Mq. In order to 
account for the total mass of hot gas (~ 4 x 10^ M©; §3), the mass loading over the nebula must 
be substantial. 

The dynamic mixing is evident, right within the central cavity around R136. HST WFPC2 
images of optical emission lines clearly demonstrate the presence of evaporative flows from various 
ionization fronts into the interior of the cavity (Scowen et al. 1998). Such flows are caused by the 
lack of a pressure confinement of HII gas after hot gas has largely escaped from the GMC through 
outflows. Both the interaction between stellar winds and evaporative flows and the interface 
between the outflows of hot and HII gases are highly unstable. The dynamic mixing can naturally 
occur. 

Supernova explosions may also play a role in the formation of the 30 Dor nebula. To be 
energetically important, however, the number of explosions needs to be several tens, at least, 
if each releases typically a few times 10^° ergs of mechanical energy. Supernova blastwaves can 
sporadically shock large amounts of evaporated HII gas in the core. The relatively high pressure 
there tends to drive the heated gas out into the halo of the nebula, similar to the outflows of 
mass-loaded stellar wind materials. 

This study of 30 Dor leads us to conclude that outflows of GHRs from GMCs can naturally 
convert large amounts of molecular gas into hot gas, via photon-evaporation, dynamic mixing, 
and occasionally shock-heating. Such mass loading can greatly enhance the radiative cooling of 
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the hot gas, therefore the X-ray luminosities of the GHRs. The hot gas properties of supergiant 
bubbles and galactic chimneys can also be affected, as they are all evolved from GHRs and have 
most likely experienced the outflow phase. 

Acknowledgments — We thank Y.-H. Chu for her comments on the manuscript. This research 
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Table 1. Results of Spectral Analysis" 


Parameters 


Values 


Plioton Index 


2 (fixed) 


norm 


3.9 X 10-4 


Low temperature 


0.150(0.145-0.159) keV 


norm 


0.94 


Higli temperature 


0.81(0.73-0.89) keV 


norm 


9.7 X 10-3 


Fe L-shell 


0.71(0.62-0.80) keV 


Gaussian Widtli 


< 9.5 X 10-2 keV 


Ne Line 


1.03(0.99-1.06) keV 


Gaussian Widtli 


< 3.1 X 10-2 keV 


Si Line 


1.35(1.34-1.36) keV 


Gaussian Width 


< 2.6 X 10-2 keV 


Column density 


1.43(1.35 - 1.48) X 10^2 cm-2 



'^The spectral normalizations are in the standard 
XSPEC units: photons keV-^ cm-2 g-i ^ j^^y 
for the power law and [10~^'^/{4:TrD'^)]J nenndV for 
the two thermal plasma components, where D is 
the distance to 30 Dor, rig and ur are the electron 
and hydrogen densities (all in cgs units). The metal 
abundances are assumed to be 30% of the solar values 
for both X-ray-emitting and absorbing gases. The 
uncertainty intervals are all at the 90% confidence. 
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Fig. 1.— RHRI X-ray image of 30 Dor. The contours are at 0.6, 1.2, 1.9, 3.0, 4.5, 6.4, 9.1, 15, 
24, 39, 63, 100, 160, 256, 407, and 647 x 10~^ counts arcmin"^, above a local background of 
~ 1 X 10-3. 



Fig. 2. — The central region of 30 Dor in X-ray and Ha. The RHRI contours start at 
3 X 10~3 counts s"^ arcmin"^. The rest is the same as Fig. 1. 



Fig. 3. — Comparison of 30 Dor in Ha (green), UV (blue), and X-ray (red). 



Fig. 4. — ASCA SIS spectrum of 30 Dor. Individual spectral components are plotted separately 
(Table 1): the power law (pink), the low and high temperature thermal components (blue and red), 
and the three Gaussian emission lines (green). 
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